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The electronic structure of YBa2Cu306+a: and La2-iSra:Cu04 for various values of x has been 
investigated using hard x-ray photoemission spectroscopy. The experimental results establish that 
the cleaving of YBa2Cu306+a; compounds occurs predominantly in the BaCuOa complex leading 
to charged surfaces at higher x and to uncharged surfaces at lower x values. The bulk component 
of the core level spectra exhibits a shift in binding energy as a function of x, from which a shift of 
the chemical potential as a function of hole concentration in the Cu02 layers could be derived. The 
doping dependence of the chemical potential across the transition from a Mott-Hubbard insulator 
to a Fermi-liquid-like metal is very different in these two series of compounds. In agreement with 
previous studies in the literature the chemical potential shift in La2-i:Sra;Cu04 is close to zero for 
small hole concentrations. In YBa2Cu306+i:, similar to all other doped cuprates studied so far, a 
strong shift of the chemical potential at low hole doping is detected. However, the results for the 
inverse charge susceptibility at small x shows a large variation between different doped cuprates. 
The results are discussed in view of various theoretical models. None of these models turns out to 
be satisfactory. 

PACS numbers: 71.10.Hf, 71.30.-f-h, 73.20.-r, 74.25. Jb, 74.72.-h, 79.60.-i 



I. INTRODUCTION 



High-Tc superconductivity in cupratesi arises from 
doping of Mott-Hubbard insulators.'^ The evolution of 
the electronic structure as a function of doping concentra- 
tion and the connected metal insulator transition is still 
not well understood and may be the central problem of 
high- Tc superconductivity. Moreover the understanding 
of the metal-insulator transition of doped Mott-Hubbard 
insulators as a material class is one of the outstanding 
challenges occupying contemporary solid state science. 
Since in these correlated systems, and in particular in 
the doped cuprates, there exists a complex interplay be- 
tween lattice, charge, and spin degrees of freedom, vari- 
ous phases between the antiferromagnetic Mott-Hubbard 
insulator (or more precisely the charge-transfer insula- 
tor) and the paramagnetic Fermi-liquid-like state have 
been the subject of discussion. The prominent state of 
matter in this context is the superconducting phase. In 
the underdoped region the pseudogap phase is also de- 



tected, which has been interpreted in terms of preformed 
pairS)^ a circulating current phasej^i^i^ii or exotic frac- 
tionalized states with topological order i^i^i^° In the un- 
derdoped region in many cuprates fluctuating or static 
'stripe' phases play an important role. These phases 
are essentially incommensurate unidirectional spin and 
charge density waves. In addition, stripe-like 'nematic' 
phases have been also predicted^^ and finally observed in 
cuprateSfi^ii^ 

An interesting question which is related to the changes 
of the electronic structure of a doped Mott-Hubbard sys- 
tem near the metal-insulator transition is the evolution 
of the chemical potential as a function of electron den- 
sity n. The rate of the chemical potential shift Sfi/Sn 
is equal to the inverse charge susceptibility Xc^ ^^'^ 
connected to the static compressibility k of an electron 
liquid by the relation k — {l/n^){6n/6iJ.). Various sce- 
narios have been discussed in relation to the question of 
the doping dependence of the chemical potential shift in 
transition metal compounds. In a simple doped semicon- 
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ductor model the chemical potential should jump from 
the first affinity states to the first ionization states when 
going from n-type doped to p-type doped systems . ^^'^^ In 
the cuprates, the former states correspond to the bottom 
of the upper Hubbard band while the latter are formed 
by the Zhang-Rice singlet states^^ (top of the valence 
band) . The energy difference between these states corre- 
sponds to the charge transfer gap which is about 2 eV. 
Thus a jump of about 2 eV in would be expected when 
going from n to p-type doped systems or of about 1 eV 
when going from the undoped system to the heavily p- 
type doped compound. An alternative scenario is that 
no shift in /z should occur as a function of doping con- 
centration due to a pinning of the chemical potential at 
impurity induced states which are in the gap.^^ A further 
model for a pinning of the chemical potential is based on 
a 'microscopic' phase separation into a metallic and an 
insulating phase, e.g. stripe or nematic phases.— 

Various concrete predictions for the chemical poten- 
tial shift as a function of the electron densit y n have 
been presented in the literature (see e.g. Ref. l2Cll) . For 
an isotropic Fermi liquid the relation for the rate of the 
chemical potential shift 

Sn ^m*' Nbifi) ^ ' 

has been given,— where m*, mf,, and Ni,{ii), are the ef- 
fective mass, the bare band mass and the bare density 
of states (DOS) at the chemical potential /i, respectively. 

is a Landau parameter which represents the isotropic 
spin-symmetric part of the electron-electron interaction. 
In an uncorrelated system, Sfi/Sn would be just given by 
the inverse bare density of state 1/Nt,{^). In correlated 
systems the effective mass would strongly increase ap- 
proaching the half-filled case. Thus a very small rate of 
chemical potential shift or a very large charge susceptibil- 
ity {S^/Sn)~^ would be expected in this approximation. 
A suppression of the chemical potential shift has been 
also predicted by numerical studies in the two dimen- 
sional Hubbard model^^-^^ for Ub close to 1, where nt is 
the number of electrons in the band. Here, from Monte 
Carlo studies a relation Afi cx —{ni, — 1)^ has been ar- 
rived at. In addition, a slave boson mean- field analysis 
of the three-band model yielded a constant chemical po- 
tential up to a hole doping p — 1/8 due to a stripe-like 
phase separation»22, Finally, based on calculations in the 
t — t' — t" — J model (which includes the long range hop- 
ping integrals t' and i"), it has been predicted^ that 
6fi/6n should be small for t'=0 and should increase with 
increasing t' . 

Experimentally, the chemical potential shift can be de- 
duced from core-level shifts in x-ray induced photoemis- 
sion spectra (XPS). The reason for this is that the bind- 
ing energy of the core level states is measured relative to 
the chemical potential /j,. More explicitly, the shift in the 
binding energy of the core levels, Ae, depends primarily 
on four terms and can be expressed as, 

Ae = -A^ + AVm + KAQ + AEr (2) 



where Afi, AVm, AQ, K, and AEr represent the change 
in chemical potential, the change in the Madelung po- 
tential, the change in the valency, a constant, and 
the change in relaxation energy, respectively.^^ It has 
turned out that in the cuprates and other transition 
metal compounds, for core levels from atoms which 
do not change their valency upon doping, the terms 
other than Afi are small and therefore the chemical po- 
tential shifts can directly derived from the core level 
shifts themselves. This core level based method was 
first applied to the system Bi2Sr2Cai_a;Y^Cu208-i-5 and 
a large total shift of the chemical potential of « 0.8 
eV between undoped and overdoped samples was de- 
duced, together with an inverse charge susceptibility 
I^Ip^o ~7 eV/holc for low hole concentration pJ^ For 
hole concentrations, p>0.1 this value was seen to de- 
crease to « 1 eV/holeJ^ Similar work has been per- 
formed on various other cuprates such as La2-2;Sr2;Cu04 
(Refs.[il|M,^,28.29,3Q), Nd2^^Ce^Cu04 (Refs. d6ij31j), 
Bi2Sr2Cai_,R,Cu208+5 (R - Pr, Er) (Refs. 
Ca2-.Na,Cu02Cl2 (RefM), and BisSrs-.La.CuOg+s 
(Ref. [IB)- Since the electronic structure close to 
the Fermi level in all these compounds is determined 
by doped Cu02 layers, a universal behavior is ex- 
pected. However, this is not observed for the p-type 
doped cuprates, since La2-a;Sra;Cu04 shows a negligible 
I^Ip^Oj but all other systems show a large |^|p^o ~1.8 
eV/holC)^ This difference has often been explained in 
terms of a 'microscopic' stripe-like phase separationji^ 
and by the fact that La2-£!;Sra;Cu04 or similar com- 
pounds in which Sr is replaced by Ba, or La is replaced by 
Nd or Eu, are much more susceptible to static stripe for- 
mation than the other cuprates . Alternative expla- 
nations of this puzzle have also been discussed in the lit- 
erature, e.g. the influence of the opening of a pseudo gap 
on the doping dependence of the chemical potential. Fi- 
nally, based on the above mentioned calculations in the 
t~t' ~t" ~ J model it was argued that for La2-a:Srj,Cu04, 
l^lp^o should be small due to its small value of t' , while 
t' for the other cuprate superconductors is larger thus 
yielding larger |^|p^o values. 

Introducing methods for the determination of chem- 
ical potential shifts, one should also mention that not 
only core-levels but also the Cu 3d^ satellite linc^® and 
the lower Hubbard band^ have been used as an intrin- 
sic reference system for the determination of the doping 
dependence of ^. 

In this article we apply the core level binding en- 
ergy method to study the chemical potential shift as a 
function of doping concentration in a further hole-doped 
cuprate, namely YBa2Cu306+2: and compare the results 
with our own data for the chemical potential shift of 
La2-a;Sra;Cu04. YBa2Cu306-i-£c is an interesting system 
since it has not only two dimensional Cu02 planes which 
are doped but also one dimensional Cu02+£c units (see 
below), which transform into metallic ribbons upon in- 
creasing the O concentration from 6 to 7. Furthermore, 
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FIG. 1: (color online) crystal structure of YBa2Cu307 (or- 
thorhombic; space group Pmmm). 0(1) sites are empty in 
the crystal structure of YBa2Cu306 leading to a tetragonal 
structure. 

in this context one should mention that YBa2Cu306+a; 
is the only compound in which nematic stripe-like order 
has been detected in the bulkJ^ 

YBa2Cu307 forms in orthorhombic structure as shown 
in Fig. 1. As evident in the figure, Y layers are sand- 
wiched by two Cu02 planes (a, 6-plane) whereby the lat- 
ter are believed to be responsible for superconductivity. 
The oxygens along the a and b axes in the Cu02 planes 
are denoted by 0(2) and 0(3), respectively and are very 
similar. The copper atoms in the Cu02 planes are de- 
noted by Cu(2) in the figure. The apical oxygens that — 
together with the Ba atoms form BaO layers — are de- 
noted by 0(4) and the oxygens in the CuO chains along 
the b directions are denoted by 0(1). The 0(1), 0(4) 
atoms and the Cu(l) atoms together form a one dimen- 
sional CuOa metallic ribbon. The decrease in x from 
x=l to x=0 leads to the fragmentation of the chains 
due to the absence of 0(1) atoms and the orthorhombic 
structure gradually changes near a;=0.3 to a tetragonal 
structurei^!^ Several studies have indicated ordering of 
oxygen atoms in the chain and/or chain fragmentsj ^^'^^i"^° 

The other end member, tetragonal YBa2Cu306 is an 
antiferromagnetic insulator. In this structure, the Cu(2) 



atoms are formally divalent, which corresponds to a sin- 
gle hole in the Cu 3d shell. In order to yield charge bal- 
ance, the Cu(l) atoms need to be formally monovalent, 
i.e., their Cu 3d states are completely filled. Evidence 
of such Cu^+ atoms has been observed experimentally 
in various studies, e.g. by x-ray absorption spectroscopy 

(XAS):41 

Increasing x from to 1 , i.e., filling the 0(1) sites, first 
the valency of the Cu(l) atoms is changed from monova- 
lent to divalent. This is completed near x—0.3. Further 
increase of x leads to the formation of formally trivalent 
Cu atoms in the chains and in the planes. Due to the 
large on-site Coulomb interaction U ~8 eV for two holes 
on a single Cu site and the smaller charge transfer en- 
ergy A «3 eV, the holes are — in reality — not formed 
at the Cu sites but reside on the adjacent O sites,— 
thus forming a Zhang- Rice singlet state together with the 
already existing hole on the Cu sites. After an insula- 
tor to metal transition, the transition into the supercon- 
ducting phase is subsequently observed. Interestingly, a 
plateau with superconducting transition temperature Tc 
of about 60 K for 0.6 < x < 0.8 has been observed, which 
is significantly different from other cuprates. Finally, a 
maximum Tc of about 91 K is observed for x close to 
0.94i 

As a reference system for YBa2Cu306+a: we have stud- 
ied La2-a;Sra;Cu04, which has the K2NiF4 structure, 
where the Cu02 planes are separated by two LaO layers. 
La2Cu04 is an antiferromagnetic insulator. Divalent Sr 
substitution at the trivalent La sites leads to hole dop- 
ing, thereby driving the system metallic. Higher doped 
compositions near x=0.15 exhibit superconductivity at 
about 40 K. 

Photoemission is usually a highly surface sensitive 
methodj^S, For photons having an energy of 1.5 keV ( 
a commonly used x-ray energy in laboratory studies) the 
mean free path of the photoclcctrons they create is of 
the order of 20 A. Increasing the photon energy to higher 
values reduces the surface sensitivity and leads to results 
which are more closely related to bulk properties. In 
addition, the knowledge of which lines are related to the 
surface and which are related to the bulk is of fundamen- 
tal importance for XPS studies on the chemical potential 
as a function of doping concentration. Therefore a pre- 
requisite for our studies of the chemical potential shift 
was a detailed investigation of the surface sensitivity of 
XPS spectra of YBa2Cu306+2: using various photon en- 
ergies and emission angles. 

Cleaving YBa2Cu306-i-x crystals perpendicular to the 
c axis can be performed at 3 different planes leading to 
6 different possible surface terminations. At present it 
is not absolutely clear whether the crystals are cleaving 
predominantly along one plane or whether other cleavage 
planes or even a mixture of the three possible cleavage 
planes is occurring which may lead up to 6 different ter- 
minations on the same cleavage plane. The appearance 
of specific surface layers may also change as a function 
of the doping concentration x. Previously, XPS studies 
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on the surface of YBa2Cu306+2; have been pubhshed for 
single crystalSj'^*'^^''*^ polycrystaUine samples,''^'''* and 
thin films.— In this context we also mention scanning 
tunneling microscopy (STM)^^. '^^'^'^i^^ studies on single 
crystals. More recently, information on the surface of 
cleaved YBa2Cu3 06+x single crystals have been also de- 
rived from Angle-Resolve Photoemission Spectroscopy 
(ARPES) studieSi^i^iSi58i59i60^^,63,64 ^Qg^ ^Yiese 

investigations of the surface layers of cleaved single crys- 
tals indicate that BaO and CuO layers are at the surface, 
i.e., the cleavage occurs predominately between the BaO 
layers and the one dimensional CuO chains. Some au- 
thors also discussed a mixed CuOs/BaO surface termi- 
nation!^ On the other hand in theoretical studies on the 
basis of density functional calculations (Refs. [H andls^ 
it was pointed out that breakage of the strong covalent 
Cu(l)-0(4) bond would be difficult. STM sudiea^ also 
reported Y layers at the surface which would indicate a 
cleavage between the Cu02 planes and the Y layers. Fi- 
nally, cleavage between the BaO and Cu02 layers have 
been reported.-- Recent ARPES studies have also de- 
tected overdoped surfaces on cleavage of optimally doped 
bulk crystals. Such surfaces with doping concentrations 
different from the bulk are expected when polar surfaces 
are produced after the cleavage process. 

In this article we present first a detailed XPS study 
of the electronic structure of cleaved YBa2Cu306-i-a; and 
La2_a;Sra;Cu04 single crystals. This is followed by the 
main focus of this contribution, the investigation of the 
chemical potential shifts in these systems as a function 
of their Cu02 plane doping concentration. 



II. EXPERIMENTAL: 

High-quality, single-phase YBa2Cu306-i-£c crystals were 
synthesized by the solution growth method^ or in non- 
reactive BaZr03 crucibles by a self flux method.— The 
crystals were examined by various other techniques such 
as inelastic neutron scattering ) muon spin rotation 
(/iSR)^ and ARPES ,^1^1^ In Table I we present the 
doping concentration x, the hole doping p in the Cu02 
layers, the superconducting transition temperatures Tc, 
and the widths of the superconducting transition tem- 
perature ATc. The non-trivial hole doping p for Cu02 
layers was derived from a recent c-axis lattice parame- 
ter study of YBa2Cu306+2; crystals*§S These values are 
in qualitative agreement with values derived from XAS 
measurements on the 01s level on untwinned single crys- 
tals.*'^ The hole doping concentration for the x=0.45 and 
0.6 samples have been determined independently to be 
p=0.085 and 0.12, respectivelyJ^ii^ La2_a;Sr;j,Cu04 crys- 
tals were grown with the traveling solvent zone method.^" 
The doping concentration of the crystals together with Tc 
and ATc are also listed in Table I. 

Despite the use of hard x-ray excitation to maximize 
the bulk sensitivity of our measurements, it is still appro- 
priate to generate a clean surface prior to measurement, 



TABLE I: Doping concentration x, hole doping p in the Cu02 
layers, Tc and ATc values of the YBa2Cu306+:c (YBCO) and 
La2-a;Sra;Cu04 (LSCO) single crystals used in this work 
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0.30 


0.30 








and so the single crystals were cleaved in a preparation 
chamber under a vacuum of lO^'^-lO"^ mbar at room 
temperature by means of removal of a top-post immedi- 
ately before an in-vacuo transfer to the analysis cham- 
ber (base pressure 3x10"^" mbar). Only the La2Cu04 
and Lai.7Sro.3Cu04 samples (i.e., two out of the total 
of eleven measured doping levels) proved to be uncleav- 
able by this method and were scraped in the preparation 
chamber using a diamond file prior to the measurements. 
We observed that the surface cleaning was necessary to 
obtain clean spectra even if hard x-rays were used for 
spectroscopy. 



The XPS measurements were performed at room tem- 
perature using the HIKE experimental station at the 
double crystal monochromator KMC-1 beamline^i at 
BESSY, Germany. This beamline operates on a bend- 
ing magnet and the energy can be scanned in the range 
of 1.7 to 12 keV. The HIKE end statioii^^ is equipped 
with a Gammadata Scienta R-4000 hemispherical elec- 
tron energy analyzer modified for high transmission and 
kinetic energies up to 10 keV. The samples were stud- 
ied using excitation energies of 2010 and 6030 eV with a 
total energy resolution of 0.33 and 0.18 eV, respectively. 
The measurements were taken at grazing incidence of the 
x-ray beam, while the axis of the spectrometer lens was, 
within a few degrees, normal to the sample surface. In or- 
der to vary surface sensitivity, off-normal emission spec- 
tra were also collected in some cases. Energy referencing 
was carried out with a combination of measurements of 
the Fermi cutoff and Au 4/ levels of a scraped Au foil, 
held in electrical contact with the samples. The accuracy 
(and reproducibility) of the absolute and relative energy 
scales is of the order of 50 meV. 
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III. RESULTS 
A. YBa2Cu306+^ 

In Fig. 2, we show the Y 3d core level spectra of 
YBa2Cu306+a; for the available x values. The spectra 
collected using 2010 eV photons in normal emission ge- 
ometry are shown in the left panel of the figure. All 
the spectra exhibit a two peak structure typical of 3c? 
core levels with spin-orbit splitting. The lines superim- 
posed on the experimental spectra are the result of a fit 
comprising two peaks shown by the thin blue lines below 
each spectrum. The intensity ratio of the two component 
peaks is close to 3/2 as expected from the multiplicity of 
the spin-orbit split features. The component peaks in 
YBa2Cu306.i5 are symmetric and become asymmetric 
with increasing oxygen concentration. This is the re- 
sult either of more than one electrostatic/screening en- 
vironment for the Y ions upon hole doping of the Cu02 
planes between which they are sandwiched, or via a mech- 
anism involving coupling of the core hole to a continuum 
of electron-hole excitations of the charge carriers in the 
Cu02 planesi^^ 

In addition to this increase in asymmetry, increasing 
oxygen doping also leads to a clear intensity shift of the 
spectral energy positions toward lower binding energies. 
Between a;=0 and 1 this shift amounts to Ae=0.8 eV. The 
analysis and discussion of this shift will form the focus of 
Sec. IV B, dealing with the chemical potential shifts as 
a function of doping. Readers for whom the description 
of the data from the different core levels — in terms of 
spectral form and the origin of the feature — would form 
an activation barrier, are free to jump directly to Sec. IV 
B and consider Fig. 15 which shows the core level shift 
data vs. doping. 

In order to test whether there is a significant surface 
component in the Y "id core level spectra, we have also 
measured the spectra at an emission angle of 35° with 
respect to the surface normal and show these data in 
Fig. 2 (b). If we would take the usual 'grazing' emis- 
sion angles of 70° or more, for the high photon energies 
used here, the signal would become impracticably small 
(due to the mismatch between the penetration depth of 
the x-rays and the inelastic mean free path length of the 
photoelectrons). Nevertheless, even the moderate emis- 
sion angle of 35° used here is sufficient to make the ex- 
periment more surface sensitive. It is evident from Fig. 
2(b) that the spectral line shapes and energy positions of 
the Y id features are essentially unchanged compared to 
those found for normal emission. This would suggest an 
absence of a surface-related Y id component. We have 
tested this hypothesis by measuring the Y id level for 
X = 0.45 with 6030 eV photons in the normal emission 
geometry, and plot the hard x-ray data (solid circles) to- 
gether with those recorded using hv=2{)lQ eV photons 
(open circles) in Fig. 2 (c). The spectra are identical, 
establishing that the Y 3d core level profile is devoid of 
surface-related features, and thus that the hv=2010eV 
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FIG. 2: (color online) Y 3d core level spectra of 
YBa2Cu306+a; for various O concentrations, (a) normal emis- 
sion spectra collected with h!/=2010eV. (b) analogous spec- 
tra collected at the same photon energy but at 35° emission 
with respect to surface normal, (c) normal emission spectra 
collected with 2010 eV photons (solid circles) and 6030 eV 
photons (open circles) for 2;=0.45. The red solid lines super- 
imposed on the experimental spectra are the total fit results 
and the thin blue lines represent the spin-orbit split compo- 
nents of the fit. 

data can be considered as wholly bulk-representative. 

Moving on from the Y 3d to the Ba 3d and 4d spectra, 
we show the 0-doping dependent data set in Fig. 3. Ba 
4d spectra collected at normal emission using hv=2{)lQ 
eV are shown in Fig. 3(a). Employment of high spec- 
tral resolution and good quality crystals reveals distinctly 
separated features in contrast to earlier data.— i^*^ The 
dominant lines at low binding energy clearly shift to lower 
binding energies with increasing oxygen concentration. 
The energy shift observed for the Ba core level lines is 
significantly larger than that observed in the Y 3d spec- 
tra shown in Fig. 2. For the Ba case the shift is Ae=1.2 
eV between a;=0 and 1. Finally, we mention that for the 
shallow Ba 5p core level excitations (not shown) the same 
shift has been detected as for the Ba 4d levels. 

The Ba 4d spectra exhibit distinct signatures of mul- 
tiple features in contrast to their Y 3d counterparts. In 
order to investigate the origin of the clear shoulders vis- 
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FIG. 3: (color online) (a) Ba 4d core level spectra 
YBa2Cu306+i for various x values collected at normal emis- 
sion with li!^=2010eV. (b) Ba Ad spectra from YBa2Cu3 06.45, 
collected at 35° emission with 2010 eV photons, and at normal 
emission with both 2010 and 6030 eV excitation energies, (c) 
A zoom in on the Ba 3d^/2 spin-orbit split component for vari- 
ous values of x collected at normal emission with h!^=2010eV. 



FIG. 4: (color online) (a) fits to the Ba 4d core level spec- 
tra of YBa2Cu306+a; for various x values collected at normal 
emission with 2010 eV photon energy, (b) Fits to the Ba 
4d spectra collected at normal emission with 6030 eV pho- 
tons, (c) The surface contribution extracted from the fits, 
as a function of x, whereby the shorthand notation 2K and 
6K has been used to designate the use of 2010 and 6030 eV 
photons, respectively. S.C. represents surface contribution. 



ible on the spin-orbit split doublets, we show the Ba Ad 
spectra for YBa2Cu3 06.45 collected at 35° emission with 
2010 eV photon energy, and at normal emission with 
hi^=2010 and 6030 eV in Fig. 3(b). In this order, the 
bulk sensitivity of the technique increases. It is clear that 
in addition to the sharp spin-orbit split features around 
87.7 and 90.5 eV binding energies, significant intensity 
grows gradually around 88.7 eV and 91.5 eV binding en- 
ergies upon increasing the surface sensitivity of the mea- 
surement. Consequently, these additional high binding 
energy features can be attributed to the core level spec- 
tra from Ba ions situated at the cleavage surface. This 
trend has also been observed for crystals with other O- 
doping levels studied in the same manner, spanning right 
across the higher doping range. The Ba 3^5/2 spectra - 
of which a portion is shown in Fig. 3(c) - also exhibit 
similar behavior with respect to the probing depth of the 
measurement. Figure 3(c) is also well suited to illustrate 
that the Ba core level line shape — for a given measure- 
ment geometry — also depends on the O concentration. 



changing from being essentially symmetric for x = 6.15 
to clearly two-component for x = 7.0. 

In order to gain a full understanding of the surface- 
bulk differences seen in the Ba Ad spectra, we have fitted 
the experimental data using asymmetric line shapes, the 
results of which are shown in Fig. 4. It is evident that 
one can simulate the Ba Ad spectra of YBa2Cu306.i5 
collected at 2010 eV and 6030 eV using one pair of spin- 
orbit split features. This indicates that the surface-bulk 
difference felt at the Ba site in this undoped composition 
is negligible. Further, we note that in every composi- 
tion, the change in surface sensitivity of the technique 
leads to significant change in the intensity ratio of the 
high to low binding energy features, thus clearly estab- 
lishing that the spin-orbit split feature at higher binding 
energies is indeed from the sample surface and the lower 
binding energy features are from the bulk. In such a case 
photoemission intensity can be expressed as: 

/(e) = [1 - e-"^/^]/" (e) -|- e-'^^^l\e) 
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where d, A, /*(e), and /''(e) is the depth of the surface 
related layer as shown in Fig. 4, the effective inelastic 
mean free path length, the surface spectrum, and the 
bulk spectrum, respectively. The first and second terms 
on the right hand side of the above equation represent 
the surface and bulk contributions to the photocmission 
spectral function. Various XPS studies indicate that d/X 
is about 0.5 for electrons having a kinetic energy of about 
1.5 keV (A '--^ 20A), which is also similar to that found in 
the universal curve i^^i^i^^ Since, A oc VEkin {Ekin is the 
kinetic energy of the photoelectrons) , d/X for Ba 4d pho- 
toelectrons is about 0.44 for the spectra corresponding 
to 2010 eV photons and 0.25 for the spectra correspond- 
ing to 6030 eV photons. Thus, the surface contribution 
in the Ba 4d spectra reduces from 36% to 22% with the 
increase in photon energy. In Fig. 4(c) we observe that 
the surface contribution (S.C.) gradually increases with 
increasing x and becomes 25% for x=l. This suggests 
that the cleaved surface contains other contributions in 
addition to BaO layers, however, the probability of BaO 
layers at the surface enhances with increasing x. 

The 1s core level spectra recorded for various O con- 
centrations are shown in Fig. 5. All the spectra exhibit 
three distinct features marked in the figure as A, B and C. 
With the reduction in surface sensitivity of the technique 
(i.e., on going from 35° emission with 2010 eV excitation 
energy, to normal emission for the same photon energy, 
to using 6030 eV photons at normal emission), the in- 
tensity of the feature C reduces significantly for all the 
oxygen doping levels studied. This indicates that feature 
C is associated with surface oxygens. 

In order to investigate the origin of multiple features 
in the O Is spectra in more detail, we have simulated 
the experimental data using asymmetric line shapes rep- 
resenting the features A, B and C. It is evident already 
from a simple visual inspection of the raw data that a 
minimum of three features are required to fit the whole 
spectral region for all the compositions. Thus, to avoid 
uncertainty in a parameterized fitting procedure, we have 
adopted a three peak fit guided by a least squares opti- 
mization with respect to the raw data. Two typical cases 
of fit results are shown in Fig. 6(a) and 6(b), illustrative 
for all the fits carried out. The doping dependent change 
in the areas of the various features is shown in Fig. 6(c). 
Firstly, we observe that the intensity of feature A re- 
mains almost unchanged in all the spectra. Feature C is 
related to the surface oxygens as argued above from its 
surface sensitivity dependence, and gradually increases 
with increasing x. The intensity of feature B also grad- 
ually increases with increasing x, which indicates a role 
for the 0(1) chain oxygen atoms in the origin of this core 
level feature. 

The binding energy of the various oxygen atoms 
in YBa2Cu3 07.o have been studied in several 
calculation a^^'^^'^^ and using XAS.**^ From these 
earlier investigations it emerges that 0(4) (see Fig. 1) 
has the lowest binding energy. The binding energies for 
the couple 0(2) and 0(3) and lastly 0(1) are about 0.7 
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FIG. 5: (color online) O Is core level spectra of YBa2Cu306+a; 
for various x values, collected using the photon energies (and 
emission angles) indicated, whereby NE stands for normal 
emission. 

eV and 0.3 eV, respectively relative to 0(4). Therefore, 
feature A in Figs. 5 and 6 which appears at the lowest 
binding energy can be attributed to 0(4) from the bulk 
of the crystal. No distinct signatures of the planar 
oxygen ions 0(2) and 0(3) and the chain oxygen 0(1) 
could be resolved in the spectra and thus, all three of 
these O sites contribute to the intensity of peak B. This 
is perfectly reasonable with regard to the data shown 
in Fig. 6 (c), in which the intensity of feature B grows 
gradually with x, as can be expected if it is partially 
0(l)-related. 

In the final part of this subsection, we present the Cu 
2p core level spectra. In the case of YBa2Cu306+a:, a 
simple attribution of each spectral feature to a particu- 
lar atomic site is severely complicated by a number of 
factors. Firstly, the existence of two types of Cu site 
— plane and chain — which take on differing valencies 
and have differing dimensionalities as a function of O 
content. Secondly, it is well known that the strong fi- 
nal state effects visible in the Cu 2p spectra of (nearly) 
divalent cuprate networks are both dependent on the 
doping level and dimension of the Cu-0 structure in 
question,^i2ii22iS 

The Cu 2p spectra collected using 2010 eV photons 
at normal emission are shown in Fig. 7. The data ex- 
hibit four spectral features arising from the spin-orbit 
split 2p3/2 and 2pi/2 levels and their satellites. In Fig. 
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FIG. 6: (color online) O Is core level spectra of (a) 
YBa2Cu306.45 and (b) YBa2Cu307.o collected under condi- 
tions giving differing surface sensitivity. For the panels con- 
taining spectra; from top to bottom the bulk sensitivity is 
increasing. The red (blue) lines represent the simulated total 
(partial) spectral functions, (c) The areas (i.e. intensities) of 
the features A, B and C as a function of x. 

7, the data have been normalized to the area under the 
2p3/2 satellite. Since the spectral line shape of Cu 2pi/2 
is very similar to the line shape of Cu 2p3/2, we concen- 
trate on the Cu 2^3/2 features in the following as shown 
in Fig. 8. In divalent cuprates, the main peak around 
933 eV binding energy represents the well-screened final 
state |2p3c?^*'L), whereby the intrinsic hole in the 3ci-shell 
has been pushed away to the ligand atoms (in this case 
O) due to the positive core-hole potentiali^l The satellite 
feature, appearing in the binding energy range of 938- 
947 eV corresponds to |2p3(i^), and is often referred to 
as the poorly-screened final state. The fact that, for the 
main line, the core-hole screening takes place via charge 
transfer to the ligands (and also thereby to neighboring 
Cu04-plaquettes^) means that an analysis of this kind of 
data offers the opportunity not only to determine param- 
eter values for the model Hamiltonians for the systems 
in question, but also to study the dynamics of charge 
transfer in cuprate networks 

Fig. 8(a) shows an overlay of a zoomed region of 



FIG. 7: (color online) Cu 2p spectra of YBa2Cu306+a; for 
various x values collected at normal emission using 2010 eV 
photon energy. 

the spectra showing the Cu 2p3/2 features for the all 
doping levels studied. Starting with YBa2Cu306.i5, as 
outlined above, Cu(l) is monovalent and Cu(2) is di- 
valent. XPS measurements on the monovalent CU2O 
compound^'' show only the well-screened peak and no 
satellite since the 3d shell is already filled in the initial 
state. Thus the satellite peak for x—0 is only due to the 
Cu(2) atoms in the planes while the well-screened peak is 
caused by a sharp Cu(l) peak and a broader Cu(2) peak 
from the planes. The broadening of the Cu(2) peak is re- 
lated to various screening channels, more localized ones 
which come at higher binding energies and more delocal- 
ized ones which appear at lower binding energy . ^^i^'^i^^ 
Increasing x leads to an oxidation of the Cu(l) atoms 
thus increasing the number of divalent Cu at the expense 
of monovalent Cu. Thus the sharp feature within the 
well-screened peak is reduced and the satellite intensity 
increases due to divalent Cu in the chains. With increas- 
ing X, similar to the Y and Ba core levels discussed above, 
the binding energy of the satellite moves to lower binding 
energies. In the context of the chemical potential shift, 
this will be discussed in Sec. IV B. The shift of the well- 
screened peak as a function of x is essentially impossible 
to evaluate, since various screening channels lead to a 
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FIG. 8: (color online) (a) Cu spectra from 

YBa2Cu306+i: for various x values recorded with 2010 eV 
photons, (b) Analogous data for the doping levels x — 6.0, 
6.45 and 7.0 collected with hi/=2010 eV and 6030 eV. 

complicated structure and the screening conditions also 
probably change for different x values. 
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FIG. 9: (color online) (a) La 3d^/2 feature and its satellite 
of La2-a;Sra;Cu04 for various x values are shown for 2010 
eV (open circles) and 6030 eV (solid circles) photon energies. 
The solid line represent the 6030 eV spectra shifted by about 
0.2 eV towards higher binding energies, (b) La 3^5/2 spectra 
using 2010 eV photons for different values of x. 



B. La2-i^Sra:Cu04 

In Fig. 9, we show the La 8(^5/2 spectra from the dif- 
ferent doping levels collected using 2010 eV and 6030 eV 
photon energies. Each spectra exhibits two features cor- 
responding to the well-screened main peak at about 833 
eV and poorly screened satellite peak at about 837.5 eV 
binding energies. The line shape of the features does 
not change with the change in incident photon energy. 
In order to verify this, we have superimposed the high 
photon energy spectra (solid line in Fig. 9(a)) by shifting 
the spectra by about 0.2 eV toward higher binding ener- 
gies. Evidently the two spectra are very similar. Since 
the bulk sensitivity of the 6030 eV spectra is significantly 
higher than the 2010 eV spectra, the above observation 
suggests that the bulk La 3d features have lower binding 
energies compared to the surface La 3ci spectra. This is 
very similar to the observations in the Ba core level spec- 
tra in YBa2Cu306+2; samples. This also indicates that 



the sample surface contains significant amount of La. 

In Fig. 9(b), we compare the La 3^5/2 features cor- 
responding to different doping levels. The vertical lines 
indicate a guide to the energy shift. The spectra at dif- 
ferent doping levels exhibit a small shift in peak positions 
indicating finite chemical potential shifts as observed in 
earlier studies. In order to verify this, for a La core level 
with lower intrinsic life time broadening we show the La 
4d spectra in Fig. 10. Each spectrum exhibits multiple 
features. To delineate the origin of various features, we 
have fit the experimental spectra using a set of doublets, 
where each doublet represents a main peak and corre- 
sponding satellite feature as seen in the La 3d spectra. 
All the spectra could be simulated remarkably well by 
two doublets with an intensity ratio of 2:3, which is sim- 
ilar to that expected for the spin-orbit split Ad features. 
Interestingly, the peak positions estimated in this way re- 
veal an energy shift as a function of doping concentration 
consistent with the La 3d spectra. 
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FIG. 10: (color online) La 4d spectra of La2-iSriCu04 for 
various values of x. The open and closed circles represent 
spectra collected using 2010 eV and 6030 eV photon energies. 
The solid line superimposed on the experimental data repre- 
sent the fit comprising a pair of spin-orbit split doublets, each 
of which is shown by a green/gray shaded surface 

The O Is spectra collected for various compositions 
and different bulk sensitivities are shown in Fig. 11. Each 
spectrum exhibits a two peak structure as demonstrated 
by fits displayed in Fig. 12. The most intense feature 
appears at about 528.5 eV binding energy. A change in 
X leads to a small shift in the peak positions as observed 
in the La Sd and 4d signals. The feature at about 531 
eV binding energy, often attributed to an impurity fea- 
ture, reduces significantly with increasing bulk sensitivity 
of the technique. However, this feature does not vanish 
even if the surface is well prepared by vacuum cleaving 
of high quality single crystals and hard x-ray excitation 
is used. Thus, this again questions its relation to impu- 
rities, raising the possibility that it is an intrinsic feature 
due to surface oxygens as also observed in YBa2Cu306+i 
samples [see Fig. 12(c)]. 

In Fig. 13, we show the Cu 2p3/2 signals corresponding 
to 2010 eV incident photon energies for different values 
of X. As described above for the case of YBa2Cu306+a; 
the satellite feature corresponds to the non-screened ex- 
citation. Interestingly, similar to the La Zd and 4d and 
O Is core level spectra a small shift in binding energy of 
the satellite feature is detected with the change in x (see 
Fig. 13). 
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FIG. 11: (color online) (a) O Is core level spectra of 
La2-a;Sra;Cu04 for various x values, corresponding to 2010 
eV and 6030 eV photon energies are shown by open and solid 
circles respectively, (b) 2010 eV spectra for different values 
of X reveal an energy shift with increasing x. 

In Fig. 14 we compare the spectra obtained with 2010 
eV and 6030 eV photon energies. The spectra are nor- 
malized to the intensity of the satellite features. In this 
case there is no energy shift visible in the satellite feature. 
In addition, the satellite features remain unaffected even 
though the technique becomes significantly bulk sensitive 
at ft,zy=6030 eV, establishing again the local character of 
this feature. 

In contrast the main peak exhibits significant redistri- 
bution of spectral weight with increasing bulk sensitivity. 
Three features can be identified in the spectra and are 
marked with A, B, and C in Fig. 14. Feature C rep- 
resents locally screened final states (whereby the ligand 
hole resides in the CUO4 plaquette containing the core 
hole) and remains almost unchanged in both the spec- 
tra. Features A and B, which are caused by more delo- 
calized screening channels are changing when changing 
the surface sensitivity in the measurements, attesting to 
different Cu environments in the surface and bulk regions 
of the sample. 
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FIG. 12; (color online) fit to the O Is features of 
La2-a;Sri:Cu04 for various x values obtained at (a) 2010 eV 
and (b) 6030 eV photon energies exhibiting two distinct fea- 
tures in all the cases, (c) The peak position and profiles of 
the O Is spectra in La2-a;Sri:Cu04 are compared to those of 
the two YBa2Cu307 end members. 

IV. DISCUSSION 

A. Surface electronic structure 

At the beginning of this discussion on the electronic 
structure of cleaved cuprate crystals (which starts by 
dealing with the case of YBa2Cu306+a;), we emphasize 
that in the present XPS experiments we have no spa- 
tial resolution. The diameter of the photon beam is of 
the order of 500 /xm. Thus we obtain only averaged in- 
formation from the surface which may, on a microscopic 
scale, possess differing terminations. On the other hand, 
due to the possibility of using photon energies between 
2010 and 6030 eV and by using different emission angles, 
depth sensitive information has been obtained, which 
can be used to differentiate between the surface and 
bulk electronic structure. In previous work on cleaved 
YBa2Cu306-i-a; single crystals (see e.g. Ref. it was 
pointed out that core level excitations (except Cu 2p) 
have a low binding energy component which represents 
bulk properties and high binding components which are 



FIG. 13: (color online) Cu 2p3/2 spectra of La2~iSra:Cu04 for 
various values of x measured with a photon energy of 2010 eV. 
While the satellite position shows only a small shift, the main 
peak exhibits significant change. The data for x=0.0 and 0.2 
are superimposed with data for x=0.04 (red triangles). 

caused by surface contaminations and by intrinsic sur- 
face contributions. Unfortunately the two high binding 
energy components appear almost at the same energy. 

As pointed out in Sec. Ill A, we see core level exci- 
tations from atoms where spectra show no sensitivity to 
the changes in the probing depth of the experiment, and 
others which show a dependence. Examples for the latter 
are the Ba and O core level excitations, while the former 
type is illustrated by the core lines from Y and the poorly 
screened final state satellite feature for Cu. 

The lack of a high binding energy peak in the Y 3d 
spectra enables us to conclude that the portion of the sur- 
face terminations containing Y layers is negligible. Fur- 
thermore, we have never observed any spin-orbit doublet 
for Y shifted by 1.5 eV to higher binding energies as has 
been detected in air-exposed samples.— This indicates 
that we do not have any Y-related surface contamina- 
tion. The slight asymmetry of the lines at higher doping 
levels can be fully understood by a Doniach-Sunjic like 
screening of the core holes by the metallic Cu02 layers 
which are in close proximity to the Y site. 

For the Ba 3d and Ad core level excitations we observe 
only one spin orbit doublet at low O concentrations, with 
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FIG. 14: (color online) Cu 2p-^/2 spectra of La2-iSri,Cu04 
for various x values collected using 2010 eV (solid circles) 
and 6030 eV (open circles) photons. A distinct signature of 
three features A, B and C is visible in the spectra of the main 
line region. 



symmetric lines which we ascribe to Ba atoms in the bulk. 
At higher O concentrations we see shoulders at higher 
binding energies. The question is whether this second 
component is related to an intrinsic surface component 
or to surface contamination. In a previous study, air 
exposed samples show a strong high energy component 
shifted by 1.5 eV compared to the bulk component 4^ 
This could indicate that the high-energy component in 
the spectra may be caused by carbonate or hydroxide 
formation due to reaction with species in the residual 
vacuum. On the other hand, we see no increase of the 
high binding energy component as a function of time. 
Therefore it is unlikely that this component is related to 
contamination due to reactions with a background pres- 
sure m the IQ-i" mbar range, or during the brief inter- 
val between cleavage at ca. 10"^ mbar and transfer to 
the analysis chamber. We thus assign the high binding 
energy component to an intrinsic surface layer and thus 
conclude that Ba must be in the terminating layer. From 
our measurements we cannot obtain detailed information 
on the possibly reconstructed BaO-CuO-BaO complex at 
the crystal termination after cleavage. At present it is 
also less clear why the surface layer in the crystals with 



low O concentration does not produce the high binding 
energy component, although a connection to the popula- 
tion of the 0(1) sites in the chains is a clear possibility. 
From the present measurements it follows that the cleav- 
age plane or the reconstruction that may follow cleavage 
probably changes between low x values and high x values. 
The spectra also indicate that for low x values there are 
no charged surface layers (which would produce a high 
binding energy component in the XPS spectra) while for 
higher doping levels, a charged surface layer exists. The 
latter results are in line with the detection of overdoped 
surface layers in ARPES spectra of YBa2Cu306+2; for 
higher x values i"^'^^'^^ 

The O Is spectra are also composed of a high binding 
energy component and a low binding energy component 
with a clear fine structure in the form of a shoulder at 
low binding energy. In Sec. Ill A, the signal below 530 
eV was ascribed to bulk O ions, with 0(4) giving rise 
to the low binding energy shoulder and the other three 
sites adding up to give the strongest feature at ca. 529 
eV. At high photon energy {hv^QOiO eV) and in partic- 
ular for small x values the intensity of the high binding 
energy component (binding energy >530 eV) is rather 
small. This indicates that using similar arguments as 
those discussed in the case of the Ba core levels, this 
part of the O Is spectrum is related to an intrinsic sur- 
face contribution. 

In the Cu 2p spectra, upon changing the surface sen- 
sitivity, no changes in the satellite features could be de- 
tected. This is in line with the interpretation of this 
feature in terms of atomic like excitations. Small differ- 
ences in the main lines are observed which could reflect 
small differences in the screening channels between the 
bulk and the surface. On the other hand, we point out 
that the changes of the main line of YBa2Cu306+2: upon 
changing the photon energy are less pronounced than in 
La2-2;Sr2;Cu04 and, as will be discussed in a moment, is 
even less pronounced than in Nd2+£i;Cea;Cu04. 

In the following we discuss the spectra of 
La2_2;Sr2,Cu04. The common view is that cleavage 
occurs between the LaO layers, thus forming LaO 
terminated surfaces. Interestingly, in the La core level 
spectra taken with different surface sensitivities, no 
surface component could be detected. Only a small shift 
of about 200 meV could be detected between spectra 
taken with hv = 2010 eV and 6030 eV. This shift may 
be caused by a Madelung potential at the La sites that 
may be different for the surface layer and the bulk. 

In principle, the bulk component in the O Is spectra at 
low binding energy should be a superposition of two lines 
stemming from the Cu02 plane O atoms and the apical 
O sites. Although we are performing high resolution XPS 
experiments, the two sites could not be resolved, indicat- 
ing that the binding energy of the O Is levels of the two 
sites is almost the same. The high energy component in 
the O Is spectra in La2-2:Sr2,Cu04 shows a strong surface 
sensitivity, almost disappearing in the spectra taken with 
the photon energy ft,z/=6030 eV. This means that similar 
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to the case of YBa2Cu306+2;, the high energy component 
is surface related. In the data from the cleaved crystals 
(x=0.04 to 0.2) this component is small, while for the 
scraped crystals (a;=0 and 0.3) the high energy compo- 
nent is considerably larger. The enhanced roughness of 
the scraped surface leads to larger surface contributions 
than for the cleaved surface. Thus, the larger high bind- 
ing energy feature in the 01s spectra from the scraped 
sample surfaces re-establishes its connection to the sur- 
face electronic structure. 

As mentioned previously the satellite feature in the Cu 
2p spectra do not change when changing the surface sen- 
sitivity of the XPS measurements using different photon 
energies in keeping with its atomic-like multiplet charac- 
ter. On the other hand, in the main line the low energy 
features A and B in Fig. 14 increase when measured with 
increased bulk sensitivity. These features are caused from 
more delocalized screening channels, in which the hole on 
the (oxygen) ligands moves further away, forming Zhang- 
Rice singlets. In a first approximation these screening 
channels would come from the in-plane electronic struc- 
ture of the Cu02 layers and therefore it is difficult to 
understand why there should be a difference between a 
surface Cu02 layer and a bulk Cu02 layer. On the other 
hand such an enhancement with increasing bulk sensitiv- 
ity was also observed in earlier studies, in particular for 
Nd2+a;Ce2,Cu04 (Ref. [soh . This could indicate that an 
explanation of the well-screened feature solely in terms 
of in-plane screening channels is probably not sufhcient. 
Possibly the total three dimensional electronic structure 
has to be used in the framework of an Anderson impu- 
rity model to calculate the screening channels for the Cu 
2p core hole to obtain a satisfactory description of the 
spectra^iSiiSa 



B. Chemical potential shift vs. doping level 

Now that we have successfully identified lines from core 
level excitations in the XPS spectra which are definitely 
related to the bulk properties and which are not influ- 
enced by surface effects, we can evaluate the bulk core 
level binding energies as a function of doping. Since the 
core level binding energies are measured relative to the 
Fermi level, we can extract the chemical potential as a 
function of hole concentration as was described in the 
Introduction. In the following, we discuss our results 
on the chemical potential shifts and compare them with 
theoretical calculations and previous experimental results 
published in the literature. 

Without further ado, in Fig. 15 we show a compilation 
of all of the relevant data for our single-crystal-only, bulk 
sensitive study of YBa2Cu3 06-(-a; and La2-a;Sra;Cu04. 
The top panel shows the shift in the binding energy of 
the core levels indicated relative to that of x = vs. x in 
YBa2Cu306-i-a;- The bottom panel shows the binding en- 
ergy shifts vs. the hole concentration p per Cu ion in the 
Cu02 planes. In Fig. 15(b) we also compare these data 
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FIG. 15: (color online) binding energy shift of the various 
bulk-related core levels lines as a function (a) of x and (b) of 
hole concentration in the Cu02 layers in YBa2Cu306+a; sam- 
ples. The analogous energy shifts for La2-iSri:Cu04 samples 
are also shown in (b). The solid lines represent a guide to the 
eye. 



with core level shifts derived in La2-a:Sra;Cu04 for vari- 
ous X — p values. For both systems the binding energy 
shift values are given relative to p = 0. 

In La2-a;Sra;Cu04 the core level excitations from O, 
La, and Cu atoms show no shift up to and including 
p=0.15. For higher p values a negative shift is detected 
(i.e., the binding energy of the core levels, referred to the 
Fermi energy is decreasing) which reaches a value of 0.2 
eV for p=0.3. Our results are in perfect agreement with 
the data published in Ref. [2^. The fact that the shifts 
Ae are zero for p <0.15 may be interpreted by an acci- 
dental cancellation of the four terms in Eq. (2) in this 
Sr concentration range. However, we believe that such a 
cancellation for three different atoms in such a large con- 
centration range is rather unlikely. Thus we assume that 
in Eq. (2) the Madelung term, the on-site valency term, 
and the relaxation term are negligible for the particular 
core level features under consideration, and that the ob- 
served shift is determined by the chemical potential shift 
A/x. We note here that, contrary to previous investiga- 
tions, we have determined the energy shift for the Cu 2p 
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lines not from the main line but from the satellite and 
thus avoid 'contamination' of the core level shift values 
with chemical shift terms related to the changing cop- 
per valence. For the O atoms the situation is less clear 
since the holes which are formed upon doping with Sr arc 
situated mainly on the O sites, as a result of the forma- 
tion of Zhang- Rice singlets. Thus, the 2p count on the 
O sites will be changed, and a chemical shift of the O Is 
core level should be expected. Apparently, the constant 
K in Eq. (2) is too small to enable the detection of this 
shift. The calculated Madelung terms in a point charge 
model are of the order of 2 eV for doping concentrations 
up to x=0.2 (see Ref. 113) • Opposite shifts are derived 
from these calculations for Cu vs. La and O core level 
excitations. On the other hand, these calculated values 
are probably strongly reduced by a high background di- 
electric constant. This assumption is supported by the 
experimental fact for the core levels of the three types 
of atoms, the same shift (or no shift) is observed, inde- 
pendent of the distance to the dopant atoms (here Sr). 
Turning to the last term in Eq. (2), it is plausible that 
the relaxation energies do not significantly depend on 
the doping concentration because the observed core level 
shifts for atoms in metallic Cu02 layers and insulating 
LaO layers are shown to be the same. We also empha- 
size that the three different types of atoms are bonded 
quite differently. For La and O in the LaO planes a pre- 
dominantly ionic bonding is a reasonable starting point, 
whereas Cu and O in the Cu02 planes are expected to 
show a strong covalent bonding component. Thus the re- 
sults for La2-£i;Sra;Cu04 would indicate that the chemical 
potential is pinned by impurity levels within the charge- 
transfer gap when going from the undoped insulating sys- 
tem to the doped metallic system. More precisely, the in- 
verse charge susceptibility 6fj,/6p is zero for < p < 0.15. 

For YBa2Cu306+2; the results for the chemical 
potential shift are completely different to those of 
La2_3;Sr2;Cu04 [see Fig. 15(b)]. Interpreting the mea- 
sured core level shifts in terms of a chemical potential 
shift, we obtain from the data for p < 0.1, presented in 
Fig. 15(b), an inverse charge susceptibility |^|p^o be- 
tween 5.5 and 9.5 eV/hole depending which core level 
lines one considers. This would support the simple semi- 
conductor scenario in which upon doping the Fermi level 
moves from the center of the gap rapidly into the va- 
lence band (effective lower Hubbard band in a single 
band model). At higher hole concentration {p > 0.15) 
the shift seems to saturate. For one group of core level 
excitations the maximum shift saturates at -0.8 eV while 
for others a maximum shift of -1.2 eV is detected. The 
lower value would correspond reasonably to half of the 
gap between the valence band and the upper Hubbard 
band in the insulating system. The higher saturation 
value is slightly too big when compared to half of the 
gap. This could indicate that the values from the Ba id 
and the O ls(B) excitations are in some way less directly 
related to the (bulk) chemical potential shift. It is re- 
markable that these two excitations both stem from the 



BaO(4) plane while the other excitations predominantly 
stem from the Cu02-Y-Cu02 complex [except the small 
contribution from the 0(1) atoms in the CuO(l) chains]. 
Having said that, it is still unclear as to why the behavior 
of the core level shifts should group in this manner. A 
priori we would have expected that the shifts from the 
ionic atoms such as Ba and Y and from the Cu 2p satel- 
lite would be the same. The large shifts of the Ba core 
level excitations as a function of doping were also dis- 
cussed in Ref. H^, however, without a possible chemical 
potential shift being taken into account. Both the Ba and 
the 0(4) atoms are close to the 0(1) dopant atoms and 
possibly experience an additional negative shift due to a 
Madelung term caused by the increasing number of 0(1) 
atoms. Since the Cu(l) atoms contribute only 1/3 to the 
Cu 2p spectra this effect is possibly not detected in the 
Cu 2p3/2 spectra. In this context one should also men- 
tion that in these Cu 2p^/2 spectra, the Cu(l) atoms only 
contribute to the satellite for p larger than « 0.3 since 
for lower p values they are monovalent and therefore do 
not form a satellite feature. 

In view of the discussions above, we tentatively assign 
the chemical potential shift in this compound to values 
measured via the Cu 2p, Y 3d, and O Is (line A) core 
level excitations. We then come to a value for the inverse 
charge susceptibihty j^\p^o ~ 5.5 eV/hole. 

This non-zero value is in agreement with values de- 
rived in previous studies. In Ref. [3 an inverse charge 
susceptibility ^|p_+o ~ 7 eV/hole has been derived 
for the system Bi2Sr2Cai_2;Y2,Cu208+5 which saturates 
for p > 0.1. The total shift between p = and 
p — 0.25 amounts to A/i — 0.8 eV. On the other 
hand in Ref. 35 a 'universal' inverse charge susceptibility 
l^lp^o ~ 1.8 eV/hole has been observed for various non- 
La2-a;Sra;Cu04 cuprates [Bi2Sr2Cai_2,Rj;Cu208+5 (R — 
Pr, Er), Ca2-a;Na:ECu02Cl2, and Bi2Sr2-2;La2;Cu06+5]. 
Thus we are left with two puzzles (i) why is 
La2_a;Sra;Cu04 the only p-type doped cuprate which 
shows a negligible inverse charge susceptibility for small 
p and (ii) why are the inverse charge susceptibilities of 
the non-La2_a;Sra;Cu04 systems so different? 

To begin with we deal with the second puzzle. For the 
various cuprates, the energy bands near the Fermi level 
originating from the Cu02 planes are fairly similar.''^ In 
this way we could understand why our values \p^o and 
A/i between p=0 and p = 0.2 for YBa2Cu306+2; agree 
reasonably with the values for Bi2Sr2Cai_2;Y-i;Cu208+5 
presented in Ref. [l^. On the other hand, it is difhcult to 
understand why YBa2Cu3 06+a; has such a different in- 
verse charge susceptibility when compared with the 'uni- 
versal' curve derived in Ref. [35l . One way to explain this 
discrepancy would be taking into account the formation 
of Cu03 chains, only existing in YBa2Cu306+a:- On the 
other hand, for small x values we start with the CuO(4)2 
dumbbells, which upon doping, i.e., insertion of 0(1) 
atoms, form more and more fragmented chains which 
cause localized states, also close to the Fermi levelj^ 
Upon completely filling the chains, a Cu(l) 3dy2_^2- 
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0(1) 2pj,-0(4) 2pz band is formed which is similar to 
the Cu02 layer bands. In addition, close to the Fermi 
level another chain band occurs which is caused by a 
TT hybridization of O 2p orbitals of the 0(1) and 0(4) 
atoms. Both this band and the localized states formed 
by the fragmented chains would be more likely to cause 
pinning of the chemical potential leading to a smaller 
chemical potential shift but not an increase of the in- 
verse charge susceptibility. Therefore it is difficult to un- 
derstand why the inverse charge susceptibility detected 
in YBa2Cu306+2; in the present experiments is so much 
larger than in the other cuprates studied in Ref. Issl . In 
this context we also mention that this large chemical po- 
tential shift in YBa2Cu306+2; is not due to formation of 
metallic CuOa chains before the insulator metal transi- 
tion appears in the Cu02 planes. XAS measurements on 
untwinned single crystals clearly indicate that the num- 
ber of holes in both CuO units shows almost the same 
doping dependence."''^ This indicates, when considering 
the dimensionality of the two units, that the insulator- 
metal transition comes first in the planes and later in 
the chains, meaning that the observed chemical potential 
shift for low doping in the Cu02 planes is not determined 
by the chains but by the planes. Finally, in this context 
we mention measurements of the plasmon dispersion by 
The dispersion coefficient extracted from 
those measurements yield information on the momentum 
dependent compressibility of the electron liquid. For the 
two compounds YBa2Cu307 and Bi2Sr2CaCu208 almost 
the same dispersion coefficient has been detected which 
fits the YBa2Cu306-i-a;/Bi2Sr2Cai_2;Ya;Cu208+5 agree- 
ment for the chemical potential shift we observe, but 
raises questions why Er and Pr-doped Bi2Sr2CaCu2 08+5 
and La-doped Bi2Sr2Cu06+5 appears to be different in 
their carefully conducted XPS experiments. 

We have confirmed in this work earlier results for a 
zero inverse charge suscegtibility for La2-xSra;Cu04 for 
small X values (see Ref. |28I) . In addition, in our present 
work we have discovered a further system which shows 
a large inverse charge susceptibility, but which is even 
larger than in various other non-La2-a;Sra;Cu04 com- 
pounds. As pointed out in the Introduction various mod- 
els for the existence or non-existence of a finite inverse 
charge susceptibility have been proposed. 

Several theoretical calculations^SiSliS^ have predicted 
for a p-type doped Mott-Hubbard insulator the appear- 
ance of a quasiparticle like resonance at the top of the 
lower Hubbard band. This would give a jump of the 
chemical potential shift A/i=J7/2. In a charge-transfer 
insulator (which can be transformed into an effective 
Mott-Hubbard insulator) the shift would be half of the 
charge-transfer gap. These calculations also predict that 
the resonances move into the gap depending on the size 
of U relative to a critical value Uc- Such a shift could re- 
duce the inverse charge susceptibility as experimentally 
observed in the system La2-a;Sr2;Cu04. On the other 
hand, in view of the similarities of the crystal structure 
and the electronic structure of the Cu02 planes in the 



various cuprate high-Tc superconductors, it is very diffi- 
cult to understand why the parameters determining the 
energy position of the quasiparticle resonance should be 
so different for the various compounds. 

Another reason for a vanishing chemical potential shift 
upon doping a Mott-Hubbard insulator could be a phase 
separation into a metallic and an insulating phase. Mi- 
croscopic phase separation has been claimed to be de- 
tected in the stripe phases in which antifcrromagnetic in- 
sulating antiphase domains are separated by periodically 
spaced domain walls to which the holes segregate (see e.g. 
Ref. [H). Actually, as mentioned in the Introduction, a 
slave boson mean-field analysis of the three-band model 
resulted in a constant chemical potential up to p=l/8. 
At this doping concentration — corresponding to four 
lattice constants in the Cu02 plane — the cuprates are 
most susceptible to stripe formation. Below this doping 
concentration, the distance between the rivers of charge 
increase more and more, leaving the hole concentration 
within the microscopic metallic regions (the stripes) un- 
changed. This could lead to a vanishing chemical po- 
tential shift, an argument which was used for the expla- 
nation of the zero inverse charge susceptibility observed 
in La2-a;Srj;Cu04 [Ref. l33|. Above p=l/8 the common 
view is that in the stripe phase, the doping concentration 
within the stripes increases. In this scenario based on mi- 
croscopic phase separation not only the vanishing chemi- 
cal potential shift for La2-a;Sra;Cu04 for small x was ex- 
plained but also the reduced Sn/Sp for Ri_2,A2,Mn03 (R: 
rare earth; A: alkaline earth) around the charge ordered 
composition range.^ In addition, the fixed chemical po- 
tential in La2_2;Sr2;Ni04 for x <0.4 [Ref. [30ll97j was also 
discussed in this context. 

On the other hand, we see serious problems to explain 
the vanishing chemical potential shift in La2-a;Srj;Cu04 
for small x on the basis of a microscopic phase separation, 
as, in this compound only fluctuating stripes have been 
observed. Only in the systems La2-a;Sra;Cu04 in which 
part of La is replaced by Nd or Eu, or Sr is replaced by 
Ba, have static stripes been detected. Furthermore, for 
x=0.15 we are definitely in a range with strongly inter- 
acting stripes which is a region far from that discussed 
in terms of microscopic phase separation. Finally, in this 
scenario it is very difficult to explain the difference be- 
tween La2_i;Sra;Cu04 and YBa2Cu306+2; as regards the 
core level shifts since in the former system only fluctu- 
ating stripes exist while in the latter system a nematic 
stripe liquid has even been detected experimentally for 
x=0.45 (Ref. [11). 

Another explanation of the anomalous inverse charge 
susceptibility in La2-3;Sra;Cu04 is based on the role of 
the next-nearest neighbor hopping integral t' in the band 
dispersion of cuprates. Theoretical work going in this 
direction was published in Ref. [2I predicting an increas- 
ing chemical potential shift with increasing t' . This sce- 
nario is related to the van Hove singularity close to the 
Fermi surface of the two dimensional electronic structure 
of the Cu02 planes. This singularity causes a high den- 
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sity of states close to the Fermi level which in turn may 
reduce the chemical potential shift. For a small t' /t {t is 
the nearest neighbor hopping integral) causing a square 
Fermi surface, the chemical potential lies in the singular- 
ity and doping causes only a small shift. Larger absolute 
t' /t values shift the singularity below the Fermi level thus 
increasing the chemical potential shifts caused by small 
doping concentrations. On the other hand, a compila- 
tion of t' /t values based on band structure calculations^^ 
shows that t' /t for La2Cu04 is considerably larger than 
that of Ca2Cu02Cl2. This would give a smaller chemi- 
cal potential shift for the latter when compared with the 
former in clear contradiction with the experimental find- 
ings. Moreover, the band structure calculations and also 
ARPES data indicate that in all the cuprates t' /t is finite, 
shifting the van Hove singularity well below the Fermi 
level. Thus a doping concentration near 0.2 is needed 
to shift the chemical potential into the sigularity which 
could reduce the chemical potential shift. In a recent 
paper therefore, a doping dependent t' /t has been in- 
troduced to explain the small |^|p^o in La2-a:Sra;Cu04. 

For completeness we also mention discussions of the 
influence of a pseudogap on the chemical potential shift 
in doped cuprates. ■2^'^ Since the observed pseudogaps 
in the underdoped cuprates are of the order of 60 meV 
they have a negligible influence on the chemical potential 
shifts detected in this work. Furthermore the existence 
of a pseudogap in underdoped cuprates is a generic phe- 
nomenon and cannot therefore explain the differences in 
the chemical potential shifts between different cuprates, 
which were discussed above. 



V. CONCLUSIONS 

We have studied the core level spectra of 
La2-2:Sr2:Cu04 and YBa2Cu306+2; samples for dif- 
ferent values of x using high resolution, hard x-ray 
photoemission spectroscopy with variable photon ener- 
gies up to 6030 cV. In those spectra we could separate 
core level features derived only from the bulk and 
features which come from the surface. The latter seem 
predominantly caused by intrinsic surface contributions, 
i.e., contributions which are related to the cleaving of 
the surface which leads to reconstructions and charge 
redistributions at the surface, and only a small part 
seems due to surface contamination. The experimental 
results establish that the cleavage in YBa2Cu306+2: 
crystals for x >0.45 occurs predominantly in the 
BaO(l)Cu(l)0(4)2 layers which result in charged ter- 
minations. For small x values the cleavage planes are 
found to be again in the Cu(l)0(4)2 layers, thus leading 
in this case to a negligible contribution from charged 
surface terminations. 

Evaluating the binding energies of the core levels 
from the unambiguously bulk-related features we could 



determine the chemical potential shifts as a function 
of the doping concentration in the Cu02 planes. For 
La2-2:Sr2:Cu04 we could confirm previous results in the 
literature indicating no shift of the chemical potential up 
to and including a doping concentration x=0.15. In view 
of a simple picture of a doped Mott-Hubbard insulator, a 
large jump of the chemical potential of the order of half 
the gap energies, i.e., of the order of 0.8 eV is expected. 
In YBa2Cu306-i-a;, on the other hand this expected large 
shift is detected when going from the undoped material 
to the optimally doped high-Tc superconductor. The 
observed concentration dependence of the chemical 
potential shift is in good agreement with that observed 
in Bi2Sr2Cai_a;Ya;Cu2 08-i-5'^ However, this shift is con- 
siderably larger than that in three other doped cuprate 
systems presented in the literature. '^^ In particular the 
changes of the chemical potential shift per hole, i.e., the 
inverse charge susceptibility for small x values differs 
considerably between different (non- La2-a;Sra;Cu04) 
p-type doped cuprates. These contradicting results are 
discussed in terms of different scenarios which could 
influence the chemical potential shift as a function of 
doping: the influence of microscopic charge separation 
(stripe formation) , particular band structure phenomena 
such as the influence of the next-nearest neighbor 
hopping integral t' , or the pseudo gap. None of these 
scenarios yield a satisfactory explanation covering the 
detected experimental results. Therefore, although we 
embarked on this systematic, bulk sensitive study of 
the doping dependence of the chemical potential fi in 
YBa2Cu306-(-x and La2-xSr2,Cu04 in order to achieve a 
generically valid picture for this fundamental property of 
these electron systems, we are obliged to close this paper 
with the call for further experimental and theoretical 
work that still seems to be necessary to understand the 
ongoing chemical potential puzzle in the cuprates. 
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